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ABSTRACT 

We show that the combination of a weak magnetic propeller and accretion disc res- 
onances can effectively halt accretion in short period cataclysmic variables for large 
fractions of their lifetimes. This may help to explain the discrepancy between the 
observed and predicted orbital period distributions of cataclysmic variables at short 
periods. Orbital resonances cause the disc to become eccentric, allowing material to 
fall back onto the donor star or out of the system. A weak magnetic field on a rapidly 
spinning primary star propels disc material outwards, allowing it to access these res- 
onances. Numerical and analytic calculations show that this state can be long lived 
(~ 10 11 yr). This is because the magnetic propeller is required only to maintain access 
to the resonances, and not to push matter out of the Roche lobe, so that the spin 
down time-scale is much longer than for a classical propeller model. 

Key words: accretion, accretion discs - stars: magnetic fields - stars: dwarf novae 



1 INTRODUCTION 

It is known that accretion disc resonances in cataclysmic 
variables ca n remove angular momentum from the outer ac- 
cretion disc llWhitehurstl l988'l . This occurs when an element 
of the disc orbits with a frequency which is resonant with 
that of the binary orbit, typically at a ratio of three to one. 
This 3:1 resonance drives the accretion disc to become eccen- 
tric and causes the superhump modulations which are the 
denning feature of the superoutbursts of cataclysmic vari- 
ables. Disc resonances provide a second source of angular 
momentum transport which is responsible for increasing the 
mass transfer r ate, driving ordinary outbursts to becom e 
superoutbursts <|Osakilll989l:lTruss. Murray fc Wvnrl200ll) . 

In known systems, superhumps normally fade soon af- 
ter outburst. This is because the tidal torques clear material 
from the outer disc until there is no material remaining at 
the resonant radius. Resonant effects cannot switch on again 
until a new outburst drives the outer disc beyond the reso- 
nant radius once more. 

In this paper we consider the possibility that access 
to orbital resonances can be maintained almost indefinitely 
by the action of a weak magnetic propeller. It is known 
that rapidly spinning and highly magnetic primary stars 
can act as magnetic propellers in some cataclysmic vari- 
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ables. Indeed, in the case of AE Aqr it is believed that 
the entire mass transfer st ream is ejected by a partic ularly 
strong magnetic propeller |Wvnn. King fc Hornelll997f) and 
the evolution of AE Aqr jschenker. Wvnn fc SpeithJl20o4) 
and WZ Sge llLasota. Kuulkers^^Charl^n999l) may well 
have been influenced by magnetic effects. Here we consider 
the possibility that a much weaker magnetic propeller might 
drive the disc outwards as far as the resonant radius, and 
so maintain a reservoir of material at this radius. The re- 
sulting continuous deposition of angular momentum by the 
resonance would drive the growth of eccentric modes. Disc 
material would be therefore be likely to extend beyond the 
primary Roche lobe. Such material would then either es- 
cape the binary or fall back to the donor star, preventing 
accretion onto the primary star in either case. 



We are motivated to study this mechanism because it 
may modify the evolution of cataclysmic variables, especially 
if this state can be long lived. This may well be the case, 
since a much weaker torque is required than in the AE Aqr 
case. The mechanism may slow or even halt accretion onto 
the white dwarf by driving mass out of the primary Roche 
lobe. At smaller mass ratios the resonant radii move further 
into the Roche lobe of the primary star, and are more likely 
to reach the disc, so that any effect is likely to influence 
short period cataclysmic variables disproportionately. This 
is interesting because existing models are unsuccessful in 
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reproducing the observed evolutionary behaviour of short 
period cataclysmic variables. 

Binary evolution models predict a minimum orbital pe- 
riod for cataclysmic variables which is reached when the 
donor star is no longer able to shrink in response to mass 
loss. M odels consistently place this minimum at around 
65 mm llKolb fc Rittall99ftlHowell. Rappaport fc Politanol 
Il997t iKolb fc Baraffelll999r e.g.). whereas the observed pe- 
riod distribution sh o ws a sharp cut off at 78m i n (e.g . 
iRitter fc Kolbl ll99St iTappert. Aueusteiin fc Mazal 120021) . 
The same models also predict that the evolution of the or- 
bital period should slow at short periods, leading to an accu- 
mulation of systems in a sharp period spike at the minimum 
period. This feature is absent from the observed distribu- 
tion, as shown in Fig. Q 

A number of explanations have been proposed for 
these discrepancies, but none are entirely satisfactory. An 
obvious solution is to hide short period systems from 
view, but standard selection effects do not act suffi- 
ciently strongly to explain the dif ferences JPattersonlll998t 
iKing. Schenker fc Hameurvl [2002). Another possibility is 
that an additional angular momentum loss mechanism could 
act to increase the minimum period jPattersorj |l998) or a 
range of angular moment um loss rate s could blunt the period 
spike iKing et alj 120021) . However, iBarker fc Kolbl (l2002h 
show that even the loss of all the angular momentum car- 
ried by the accretion stream cannot increase the minimum 
period sufficiently. They also argue that blunting the period 
spike by 'smearing' the angular momentum loss rates would 
require a great deal of fine tuning to be successful. 

In this paper we consider how the effects of propeller- 
induced access to orbital resonances might prevent some 
short-period CVs from exhibiting outburst behaviour, and 
cause them to become fainter. In Section [5] we discuss the 
model in more detail. In Sections [3] and we explore the 
parameter space for which the idea may be of interest us- 
ing smoothed particle hydrodynamics calculations. Finally, 
in Section we discuss how this mechanism might influence 
the evolution of short period cataclysmic variables. 



2 PROPELLER-DRIVEN RESONANT MASS 
LOSS 

2.1 Orbital Resonances 

In this section we test the idea described in Section by 
establishing under which conditions the accretion disc has 
access to the orbital resonances. These conditions are con- 
veniently described with reference to a series of important 
radii, measured from the primary star, all of which are plot- 
ted in Fig. |2| Disc elements describe orbits that are approx- 
imately Keplerian, so that the positions of the orbital reso- 
nances are readily calculated from Kepler's law. Therefore, 
where a is the binary separation, the 2:1 and 3:1 resonances 
are found respectively at 

2/3 

' (i + gr 1/3 , (i) 
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a 

If an accretion disc is to have access to a resonance then 
the resonant radius must lie within the disc. If the disc ra- 
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Figure 1. CV orbital period histogram with the cumulative fre- 
quen cy overplotted as a solid line. This plot is based on fig. 1 
from lKing et aJjfcOOSlK but uses more recent data from then on- 
line version of lRitter fc Koltl ll998l) . The plot contains data for 
all 215 CVs with periods known to be between 60 and 120 mins, 
from a total of 569 with known periods in the catalogue. 
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Figure 2. Plot of resonant, Roche lobe, circularisation, and vis- 
cous disc radii about the primary star in a close binary as a func- 
tion of mass ratio q. 



dius can be related to the mass ratio then this consideration 
can be used to calculate q ma x, the maximum mass ratio for 
which resonant access is possible. There are several ways of 
estimating the extent of the disc. Firstly we use the Roche 
lobe radius. The Roche lobe of the accreting star provides 
a hard maximum for the disc radius. The disc is, in fact 
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Figure 3. This figure illustrates the positions of the 2:1 and 3:1 resonances in binaries with mass ratios of q = 0.05 in panel a, and 
q = 0.5 in panel b. The primary Roche lobe is marked with a solid line while the secondary star, which is assumed to be Roche lobe 
filling, is coloured grey. The two dashed circles represent the resonant radii, the inner circle representing the 3:1 resonance and the outer 
line being the 2:1 resonance. The length unit is the binary separation a and the diagram is centred on the centre of mass. The primary 
star is represented by a cross. In case a, both the resonant radii are within the Roche lobe and accessible to the disc, whereas neither is 
available in case b. 



unlikely to extend this far in quiescence, but this provides 
a good estimate for the radius in outburst, when the disc 
is subject to increased vi scous spreading. The Roche lobe 
radius is approximated bv lEggletonl il983T> to be 



a 



0.49g 



-2/3 



(3) 



0.6<r 2 / 3 +ln(l + q- 

where q represents the mass ratio of the binary, q = M2 jM\ . 
Equation [3] shows that the primary Roche lobe contracts 
with increasing q, as must the maximum extent of the ac- 
cretion disc. Equating the 3:1 resonant radius and the Roch 
lobe radius, using equations Q and and solving for q 
gives a maximum mass ratio for access to the 2:1 resonance 
thus 



9max,2:l = 0.06 



(4) 



The same technique when applied to the 3 :1 resonance gives 
a maximum access ratio of g max ,3:i = 0.90. |Paczvriskil jl977f) 
computed g ma x,3:i by solving the restricted three-body prob- 
lem, finding that 



g m ax,3:l = 0.25 



(5) 



This much smaller value for qmax, 3 : 1 reflects the fact that 
the accretion disc is tidally truncated within the primary 
Roche lobe. These results are also seen graphically in Fig. [5] 
as the intersections of the resonant radii and the Roche lobe 
radius. This approach therefore provides a good estimate of 
the maximum mass ratio g max for which resonance access is 
permitted in outburst. For reference, the positions of the 3:1 
and 2:1 resonances are also marked in Fig.^lfor two typical 
CV mass ratios. 



We note that q ma x,2:i is precisely th e mass ratio at- 
tributed to WZ Sge ijSkidmore et, aljr2002h . Indeed an out- 
burst orbital hump (O OH), probably dr i ven b y this reso- 
nance, was observed by IPatterson et alJ (l2002f) before the 
onset of 3:1 superhumps during the 2001 outburst, sug- 
gesting that the disc can access this resonance during out- 
burst, as would be expected. SU UMa stars exhibit super- 
humps due to the 3:1 resonance so one would expect them 
to obey con dition 0. Indeed, of all the confirmed SU UMa 
in lRitter fc Kolbl Jl998l) . only two have q > 0.25. 

A rough estimate of the quiescent disc radius is given by 
the circularisation radius R CY i c . This is the radius of a circu- 
lar orbit associated with a specific angular momentum iden- 
tical to that with which mass crosses the first Lagrangian 
point. It is the radius at which a disc will begin to form, 
and so all discs must extend at least this far. The circular- 
isation radius may be used as an esti mate for the extent of 
a ver y cold disc and is given by (e.g. iFrank. King fc Raind 
2002) 



Rcii 



(1+9) 



(6) 



Here, the position of Rli, the first Lagrangian point, may be 
found numerically, or it may be fitted by various expressions. 
The position of Rli should not be confused with the Rl, the 
mean Roche lobe radius as defined in equation © • iKooall 
lll959l) gives the following fit for Rl,i. 



Rli = 
a 

where 



1 2 1 3 



(7) 
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3 (l + q) 



(8) 



Equating the 3:1 resonant radius and the circularisation ra- 
dius using equations and @, and solving for q, gives 
g m ax,3:i = 0.01 as the maximum mass ratio for which the 
3:1 resonance is accessible, if the disc is limited by the cir- 
cularisation radius. Such a limit however implies that there 
is no outwards viscous spreading whatever. Even quiescent 
accretion discs are however subject to some viscous spread- 
ing. 

An improved estimate for the extent of the quiescen t 
accretion disc is that proposed bv lOsaki fc Meveil {2002). 
In this approach conservation of angular momentum in a 
steady-state viscous disc, neglecting external torques, is used 
to estimate the outer radius of an accretion disc thus 



Rou 



ay 



Rciv 



(9) 



for a quiescent Shakura-Sunyaev disc. Combining this result 
with equation @, this approach gives an estimate of the 
maximum mass ratio for 3:1 resonance access in the quies- 
cent disc of 



9max,3:l 



0.29 



(10) 



We note however that access to the 2:1 resonance is still lim- 
ited by the Roche lobe estimate given by equation[I]since the 
Roche lobe is smaller than the viscous limit in this regime, 
as shown in Fig. |2] 




0.20 



Figure 5. Spin-down time-scale t as a function of mass ratio q for 
a typical CV. In the resonant case, plotted with a solid line, the 
time-scale increases asymptotically as the circularisation radius 
approaches the 3:1 resonance while r sp i n ~ 10 10 yr for most q. In 
the non-resonant, or strong propeller case, plotted with a dashed 
line, the spin-down time-scale becomes large only for q < 0.005. 



2.2 The magnetic propeller 

So far we have discussed access to orbital resonances for an 
unperturbed, quiescent disc. During outburst the disc will 
expand outwards due to increased viscous angular momen- 
tum transport, hence the detection of superhumps in super- 
outbursts. Another way to aid access to the binary reso- 
nances is by the addition of an additional source of angular 
momentum, such as a magnetic torque. A magnetic field on 
the primary star can act either to drive disc matter inwards 
or outwards; it may behave either as a source or as a sink 
of angular momentum. It is useful here to examine a simpli- 
fied case where the primary magnetic field rotates as a rigid 
body, rooted to the primary star and rotating with angular 
frequency f2* . The magnetic or magnetospheric radius 7? mag 
is the radius at which the magnetic time-scale is equal to 
the viscous time-scale; in other words it is the largest radius 
at which magnetic torques are significant in the disc. The 
co-rotation radius R co is the radius at which an element of 
disc, in describing a Keplerian orbit about the primary star, 
shares the angular frequency of the star and the magnetic 
field. At this co-rotation radius, analogous to that of a geo- 
stationary orbit, there is no drag force between the field and 
the disc, since they have the same velocity. However 7? co 
marks the boundary between the magnetic accretion and 
magnetic propeller regimes, as illustrated in Fig. 2] 

At those radii where R < R co the magnetic field rotates 
less rapidly than the orbiting disc, and so acts to brake it, 
aiding accretion. If 7? mag < R co therefore, then the disc will 
be be disrupted at i? mag by the effect of rapid magnetic 
accretion, causing the star to spin up, as is the case in an 
intermediate polar. However, if i? mag > -R co then then disc 
material, spreading inwards, first encounters a significant 



magnetic torque where the field is travelling faster than the 
disc, causing the transfer of angular momentum from the 
star to the disc, pushing the disc outwards and spinning the 
star down. Although magnetic accretion and the magnetic 
propeller both result in a 'hole' in the accretion disc, they 
are distinct mechanisms with different outcomes for the pri- 
mary spin and for the disc structure. Disc truncation by 
the magnetic propeller is implicated in the long recurrence 
time of WZ Sge, and in FU Ori young stars dMatthews et alJ 
2004). AE Aqr is an extreme case of propeller truncation. 

In outbursting CVs, such as dwarf novae, an outburst 
cycle exists and a steady state is never reached. A propeller 
may reduce the disc mass sufficiently so that such an out- 
burst is never triggered and a steady state is reached. In 
this steady state, the disc must lose mass at the same rate 
at which it is fed by the accretion stream. This can lead 
to only two possible steady-state solutions. The disc will 
spread viscously as mass is added and either the inner disc 
will reach R co and the system will become a magnetic ac- 
cretor, or the outer disc will extend beyond Rl and mass 
will be lost from the Roche lobe. In the latter case, that of 
the true propeller, a strong magnetic field is required, since 
mass must be pushed to radii beyond 7?l, and a rapid spin 
down is therefore suffered by the primary star. Such a state 
cannot therefore be long lived and such systems are there- 
fore likely to be rare, and will not be significant in the CV 
period distribution. 

2.3 Combined effect of resonances and propeller 

The model proposed in this paper uses both magnetic and 
resonant torques to enable a weak magnetic propeller to 
eject mass from the primary Roche lobe. The model re- 
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Figure 4. This figure illustrates magnetically enhanced accretion in a, and the magnetic propeller in b. The primary Roche lobe is 
marked with a solid line while the secondary star, which is assumed to be Roche lobe filling, is coloured dark grey. The accretion disc 
is shaded light grey and is bounded at its inner edge by a dashed circle representing the magnetic radius. The dotted line denotes the 
co-rotation radius. The length unit is the binary separation a and the diagram is centred on the centre of mass. The primary star is 
represented by a cross. The difference between the two cases is explained fully in the text. 



quires the magnetic propeller to push mass only as far as 
the 3:1 resonance. This implies a much weaker propeller, 
and hence a longer spin down time. Such a propeller state 
may therefore exist for much longer than one which relies 
on the magnetic torque alone. In this picture, disc elements 
that are propelled to the resonant radius will be excited 
into eccentric orbits, losing angular momentum. The newly 
eccentric orbits may carry disc material outside the Roche 
lobe at apastron, where it will be ejected from the binary 
or re-accreted by the secondary star. Both the propeller and 
the resonant torque are critical to this model. 

Propeller-driven resonant mass loss will occur if two 
conditions are satisfied. Firstly, the circularisation radius 
must also lie close to the 3:1 resonance so that a weak mag- 
netic torque can lift mass to the resonant radius. Secondly, 
the 3:1 resonance should lie within, but ideally close to, the 
Roche lobe radius so that a small eccentricity is able to drive 
mass from the Roche lobe and cause mass loss. 

The magnetic torque and hence the white dwarf spin- 
down time-scale required to drive resonant access can be 
estimated by calculating the torque required to lift mass 
from the circularisation radius to the 3:1 resonant radius, 
assuming that such mass will then be lost to the Roche lobe 
as a result of resonant torques. In the steady state the rate at 
which mass is transferred from i? c irc to i?Li is identical to the 
binary mass transfer rate M. This gives a magnetic torque, 
which in the absence of accretion or any other sources of 
angular momentum, is identical to the primary spin-down 
torque, of 



Jmag = V GMlM 



R, 



1/2 



(11) 



where G represents the constant of universal gravitation and 
Mi is the mass of the primary star. Equations © and \FQ 
may be combined to give 



Jn 



VGM7A/a( 7 1/2 -[l + g] 1/2 



(12) 



where i?Li is obtained for example from equation (|7J and 
where 7 = 0.48 for the 3:1 resonance. The spin-down time- 
scale is given by 

Jwd 27r/ w d 



Jn 



(13) 



Pspin J mag 

where J w d is the spin angular momentum of the white dwarf, 
Jwd ~ MiR^ d is the moment of inertia, R vd is the radius 
of the white dwarf and P sp in is the WD spin period. The 
spin-down time-scale is then takes the form 



1/2 



P^mMG 1 / 2 



I -"3:1 - H circ 



(14) 



The spin-down time-scale is then a weak function of a and 
Mi, but a strong function of q, rising asymptotically at 
Po3 : i = i?circ- For a CV with orbital period P or b = 80 min, 
primary spin period P sp m = 30 s, P w d = 1 X 10 9 cm, mass 
transfer rate M = 3 x 10" 11 Mq yr -1 and primary mass 
Mi = I.2M0 (selected to mtach numerical simulations) the 
spin-down time-scale is roughly constant at r sp i n ~ 10 yr 
for most q, but is much longer for q < 0.05. This is illus- 
trated in Fig. with an asymptote at q ~ 0.02. For q < 0.02 
the circularisation radius lies outside the 3:1 resonant radius 
so that only the inner disc would have access to the torque 
from this orbital resonance. This may explain the results 
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of ISimpson fc Woodl il998l) who find, in SPH simulations 
of non-magnetic systems, that an eccentric precessi n g disc 
does not form where q < 0.025. ISimpson fc Woodl (|l998h 
themselves invoke the damping effect of the 2:1 resonance 
to explain this effect. 

For comparison, in the case of a pure magnetic pro- 
peller, where resonant torques are not considered, mass must 
be lifted by the magnetic torque all the way from the cir- 
cularisation radius to the Roche lobe surface. Applying the 
same method as above this gives 



(15) 



The spin down time can be computed in the same way, and 
exhibits the same asymptotic behaviour. However, as shown 
in Fig. the spin-down time exceeds 1 x 10 10 yr for the 
non-resonant case only for unphysical mass ratios where q < 
0.005. 

In the standard picture of CV evolution, a large num- 
ber of systems are expected to accumulate around the period 
minimum, with small mass ratios. The mechanism outlined 
above effectively halts accretion onto the primary star for 
a significant fraction of the CV lifetime. During this phase 
outbursts are prevented and the disc mass is greatly reduced, 
also reducing the disc luminosity and making the system ef- 
fectively undetectable. This phase will be much longer lived 
for low q CVs so that more short period CVs would be hid- 
den from view, providing a possible explanation for the ab- 
sence of the period spike. 



3 SMOOTHED PARTICLE HYDRODYNAMICS 
CALCULATIONS 

The smoothed particle hydrodynamics (SPH) method is a 
Lagrangian particle technique for the modelling of fluid dy- 
namics. The fundamental principle of SPH is that fluid prop- 
erties such as pressure and density are carried by particles 
which represent elements of the fluid in question. The val- 
ues of these properties may then, at any point in space, 
be interpolated using a weighted mean of the values car- 
ried by nearby particles. Conservation of angular momentum 
is inherent in the SPH technique, which makes it an ideal 
method for the modelling of accretion dis cs. A review of the 
SPH method in general can be found in iMonaghaij 1( 1992') 
while the part icular code used h ere is discussed bv iMurravl 
l|l996ft and bv lTruss et"aD <2000h . The viscosity used in this 
wo rk is the linear el ement of the SPH viscosity described 
bv iMonaeharJ lll992l) . In the divergence free case, this vis- 
cosity can be shown to b e equivalent to the Navier- Stokes 
shear viscosity an d to thelShaku ra fc Sunvaevl j l97.?l) alpha 
disc viscosity (e.g. |Pongracicll 988: Matthews 2004J). In three 



dimensions this viscosity is given by 



v = — (c s H 
10 s 



ac s H 



(16) 



where £ is the SPH linear viscosity parameter, a is the 
Shakura-Sunyaev viscosity parameter, c s is the sound speed 
in the disc and H is the disc scale height. In this version of 
the code the sound speed is fixed throughout the calculation. 

Full magnetohydrodynamics have not yet been success- 
fully incorporated into SPH though efforts have been made 
in this direction (e.g. IPrice fc Monaghanl 12005). A simple 



approach to modelling the effect of the primary magnetic 
field is to apply a magnetic torque prescription to the SPH 
particles. The magnetic element of the momentum equation 
can be expressed as 



J x B = — (B-V)B-V 



B- 

2/', 



(17) 



where J is current density, B is magnetic field and p T is the 
permeability of the plasma. The first term on the right hand 
side is in the form of a tension, which vanishes if the field 
lines are not curved, and the second represents a pressure 
which can usually be neglected for slowly varying fields. 

If equation 1171 1 is transformed into cylindrical polar 
coordinates, and if it is assumed that the radial component 
of the magnetic field is small and that B 2 ~ B Z B^, then th e 
magnetic acceleration takes the form (e.g. lMatthewsll2004) 



J_ f B K B v \ 
pr c V 47r / ' 



(18) 



where p is density and r c represents the radius of curvature 
of the magnetic field lines, while B z and B v represent the 
vertical and azimuthal components of the magnetic field re- 
spectively. We make the further approximation that the ra- 
dius of curvature of the field lines i s of the order of the disc 
scale- height so that r c ~ H (e.g. IPearson. Wvnn fc Kind 
Il997l) . The ratio of vertical to azimuthal field strengths is re- 
lated to the shear between the disc and the magnetic field. If 
it is assumed that the field rotates rigidly w ith the star then 
this r atio can be expressed in the form (e.g. iLivio fc Pringlel 
1992) 



(19) 



B v [Q,\ — O*) 
B z 

where fik denotes the Keplerian angular velocity at radius 
R from the star. If Keplerian orbits are assumed, the surface 
density related to the disc scale height by E ~ pH, and the 
magnetic moment is defined by the expression p = B Z R 3 , 
then we may make the approximation 

p 2 R~ 6 (fi k 



4ttE 



P 2 R-" 
4tt£ 



|Vd - vb 

Vd 



(20) 



where Vd and v_g are the velocities of the disc and the mag- 
netic field respectively. It is an acceleration of this form 
which is added to the particles in the SPH model to rep- 
resent the effect of the primary magnetic field. The acceler- 
ation is fixed in the code by the parameter 



4tt£ 



(21) 



The effects of the binary resonance require no such effort to 
reproduce, but emerge naturally as a result of gravitational 
forces exerted on the particles by the binary. Every parti- 
cle is subject to a gravitational acceleration from each star 
in addition to the magnetic acceleration defined in equation 
il'iOt and fluid dynamic forces. In order to compute the pre- 
dicted spin evolution of the primary star, the torques due to 
magnetic and gravitational accelerations are calculated and 
recorded separately. However, at this stage the code per- 
mits neither the spin of the primary nor the orbital period 
to evolve as a result of the computed angular momentum 
transfer. 

Mass is injected at the first Lagrangian point and is 
removed from the simulation in three ways. Particles which 
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reach the surface of the white dwarf at R* — 6 x 10 cm are 
removed and count towards accreted mass. If elements of 
the disc extend beyond 1 x 10 10 cm then they removed and 
counted as mass lost to the system. If mass re-enters the 
secondary Roche lobe, then it is considered to have been 
re-accreted by the secondary star. 



4 NUMERICAL RESULTS 

A series of numerical experiments were performed to con- 
firm the viability of the mechanism explained in Section 
and to briefly explore the parameter space for which that 
mechanism is operable. Certain basic parameters of the 
system were fixed: the total mass of the system was set 
to Mtot = 1-2 Mq, because the simulations were initially 
part of a programme to model WZ Sge. The orbital pe- 
riod was Porb = 80 min and the spin period of the pri- 
mary star was a deliberately short P sp m = 30 s. Particles 
were injected into the primary Roche Lobe at a rate of 
M = 3.0 x 10" 11 M Q yr _1 although in some cases little of 
this mass found its way to the primary star. The Shakura- 
Sunyaev alpha viscosity was fixed at a — 0.01 throughout 
the calculations. In all cases the computation began with an 
empty disc and was permitted to run until a steady state 
was reached. In reality such systems would evolve, with full 
discs, from longer periods but the choice of an empty disc as 
an initial condition provides a clear illustration of the onset 
of resonant mass ejection and, since such simulations can 
be performed rapidly, allows a greater number of computa- 
tions to be performed. The particle mass was varied between 
models to enable steady state to be achieved in those compu- 
tations which led to high mass discs, without causing lower 
mass discs to be under resolved. 

The experiment was performed with mass ratios of 0.05, 
0.1, 0.2, 0.4 and 0.8 and the magnetic parameter fe m was set 
to values which, for a primary star of radius P+ ~ 1 x 10 9 cm 
and a disc surface density of E ~ 1 g cm -2 , correspond to 
very modest magnetic field strengths of B z — 5 kG, 4 kG, 
2 kG and zero at the surface of the white dwarf. It should be 
noted however that the relationship between the magnetic 
field strength and the code parameter fc m is model depen- 
dent, so that results should be treated as no more than a 
guide to the real fields required to support this mechanism. 

The results showed that both the steady-state disc mass 
and primary accretion rate were an order of magnitude lower 
in the magnetic cases than in the non-magnetic case; the disc 
mass fell from 10~ 13 to 10 _ls of the binary mass while the 
accretion rate onto the primary fell from 0.9 to 0.1 of the 
mass transfer across the LI point. Neither the magnetic field 
strength amongst the magnetic cases, nor the mass ratio had 
a clear effect on the steady-state disc mass or on the white 
dwarf accretion rate. However, these simulations were per- 
formed at mass ratios higher than those for which a strong 
dependence would be expected. 

An example of resonant mass ejection, drawn from the 
range of calculations described above, is illustrated in Fig.^l 
Disc matter is prevented from accreting by the magnetic pro- 
peller and the disc spreads outwards. When the disc grows 
sufficiently to be able to access the 3:1 resonance, eccentric 
modes are excited in the disc. The eccentricity grows until a 
proportion of the disc protrudes beyond the primary Roche 
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Figure 7. This figure shows the variation of disc mass with time 
during the onset of propeller-driven resonant mass ejection, be- 
ginning from an empty disc. Results are taken from the same SPH 
computation as for Fig. |^] Extensive mass loss commences after 
40 binary orbits and a quasi-steady state is then reached. From 
this time onwards a modulation appears, at a frequency close to 
that of the binary orbit. N ote the agreement betw een this result 
and that shown in fig. 5 of lHirose fc O saki (1990). 

lobe and can be swept up by the secondary star, or is lost to 
the system. That this eccentricity is not the result of initial 
particle velocities is clear from the fact that the disc passes 
through a circular phase and that the effect is absent in the 
non-magnetic case. 

The longer-term evolution of the disc in the same com- 
putation is illustrated in Figs. and [g] where the disc mass 
and the torque exerted upon the disc are plotted against 
time for around a hundred orbits. It is clear that once the 
disc 'catches' the 3:1 resonance after about 40 orbits, disc 
mass is lost rapidly from the system until a quasi-steady 
state is reached. This result is in good agreement with sim- 
ulations performed bv lHirose fc OsakH ijlOOOJi . Both the disc 
mass and the gravitational torque upon the disc exhibit a 
modulation on roughly the orbital period, which emerges 
at the same time as the aforementioned dramatic mass loss 
episode. This is easily understood by considering the sec- 
ondary star sweeping up disc mass once per orbit from the 
protruding part of the disc, and exerting a similarly modu- 
lated resonant torque. In fact, taking the precession of the 
disc into account, the modulation of the disc mass would 
be expected to occur on the beat frequency between the bi- 
nary orbit and the disc precession period: the superhump 
period. This indeed proves to be the case. No such mod- 
ulation is observed in the magnetic torque which remains 
roughly constant with time after the quasi-steady state has 
been reached. This is as expected, because the magnetic field 
acts predominantly on the inner edge of the accretion disc. 

By definition, the total angular momentum of the disc 
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Figure 6. These plots, taken from an SPH computation, of log surface density in g cm -1 show how propeller-driven resonant mass 
ejection might begin, starting from an empty disc. The simulation contained ~ 30, 000 particles with a sound speed of c s ~ 1 X 10 s cm s — 1 . 
The primary Roche lobe is marked with a solid line while the secondary star, which is assumed to be Roche lobe filling, is coloured 
grey. The dashed circle represents the 3:1 resonant radius. The length unit is the binary separation a and the diagram is centred on the 
centre of mass. The binary is orbiting in the anticlockwise (direct) sense. The primary star is represented by a cross. This simulation 
was performed with a mass ratio of q = 0.05, a primary mass of Mi = 1.2 Mq, a binary period of P or b = 80 min, a primary spin period 
of 

Ppin — 30 s, a mass transfer rate in the accretion stream of M — 3.0 X 10 11 Mq yr 1 and a surface magnetic field of 5 kG. In 
frame a, after a single orbit, the disc is just beginning to fill and the initial rosette structure is still visible. In frame b, after 15 orbits, 
the disc has become circularised as a result of viscous processes. However, in frame c, after 25 orbits, the resonance is beginning to drive 
an eccentricity, and by frame d, after only 50 orbits, that eccentricity is sufficient to cause extensive mass loss, with a large amount of 
mass being swept up by the secondary star. 



© 0000 RAS, MNRAS 000, 000-000 



The fate of short period cataclysmic variables 9 



2 i—i — i — i — | — i — i — i — | — i — i — i — | 

Magnetic 

^ : Tota 

-5 Gravitational 



a> 
cn 
S~ 
a 

SL 




c 

< 



-2 r i 1 

20 40 60 80 100 120 

Time (binary orbits) 

Figure 8. This figure shows the change of angular momentum 
due to magnetic and gravitational torques during the onset of 
propeller-driven resonant mass ejection, beginning from an empty 
disc. Results are taken from the same computation as for Figs. l6l 
and|7| The effect on disc angular momentum due to mass flux is 
not illustrated. 

does not exhibit long term variation after the quasi-steady 
state has been reached. The gravitational torque acts to 
counteract that applied by the primary magnetic field. The 
disc is spun up by the magnetic field and spun down again 
by the binary, so that a weak spin couple exists between 
the binary orbit and the primary spin. This drives a long 
term tendency towards a synchronous lock, but on a time- 
scale so long as to be largely irrelevant. Fig.[S]shows a small 
net torque acting to spin down the disc. However, the plot 
shows only applied torques and does not include the flux of 
angular momentum carried by particles entering and leaving 
the computation, which counteracts this trend, and allows 
a quasi-steady state to be maintained. 

An approximate spin-down time-scale for the primary 
star can be calculated from equation W'M where J mag is cal- 
culated here from numerical results, disregarding accretion 
torques. A value may be taken from Fig. [S] for example. In 
all cases where accretion was prevented in the above com- 
putations, spin-down times are found to be of the order 
yr in good agreement with analytic treatment 
in the previous section for q < 0.05. Numerical simulations 
have not yet been performed for very small q values, but 
analytic results show that the spin-down time-scale would 
be much longer. 



5 EVOLUTIONARY CONSEQUENCES 

It has been shown that accretion, and hence luminosity, may 
be greatly reduced in CVs by a weak propeller acting to- 
gether with orbital resonances. In this mechanism the mag- 



netic field of the rapidly spinning primary drives matter to 
the 3:1 orbital resonance. The action of this resonance then 
excites eccentric modes in the disc material. Disc mass is 
then lost from the primary Roche lobe and some may be re- 
accreted by the donor star. A CV may exist in this state for 
a significant fraction of its lifetime. Importantly, this non- 
accreting state can last much longer in systems close to the 
period minimum, so that a larger proportion of short period 
CVs would be expected to exist in this invisible state. This 
mechanism of propeller-driven resonant mass ejection may 
therefore go some way towards explaining the absence of the 
predicted excess of short period CVs; the period spike. 

In this picture, CVs evolving towards shorter periods 
will enter the non-accreting weak propeller state at vari- 
ous mass ratios, depending on the spin and magnetic field 
of the accreting white dwarf. This state is reached as soon 
as the 3:1 orbital radius becomes accessible to magnetically 
torqued disc matter. This phase will continue until the white 
dwarf has spun down sufficiently that accretion can resume, 
which occurs on a time-scale of r sp i n ~ 10 10 yr for most mass 
ratios but which increases greatly at very small q. Because 
the non-accreting state is longer lived in short period sys- 
tems we would expect this mechanism for hiding the period 
spike to be more efficient for those CVs entering this weak- 
propeller state later in their binary evolution, i.e. those with 
weaker primary magnetic fields. If propeller-driven resonant 
mass ejection does occur in CVs then it will greatly reduce 
the luminosity due to the accretion disc and to accretion 
onto the white dwarf. A population of low-luminosity CV- 
like objects with low luminosity is therefore predicted by this 
model, for which there are currently no known candidates. 
However, since in addition to a greatly reduced luminosity, 
outbursts would also be prevented in such systems, the non- 
detection of these objects may be due to selection effects. 
If such stars were detected however, then they would be 
expected to exhibit permanent superhumps, since the disc 
would always be eccentric and precessing. It may also be pos- 
sible to observe the re-accretion of matter by the secondary 
star in these objects. 

In normal CV evolution, mass transfer from the sec- 
ondary to the primary star, driven by an angular momen- 
tum loss mechanism such as gravitational radiation, reduces 
the mass ratio with time. In the weak propeller mechanism 
proposed here, most mass is not accreted by the primary 
star, but instead is re-accreted by the secondary so that the 
mass ratio does not evolve. For the secondary Roche lobe to 
accept this mass, either the contraction of the Roche lobe 
must be slowed or the secondary star itself must be made to 
contract more rapidly. The weak propeller mechanism lends 
itself to both of these possibilities. Matter is returned to the 
secondary with increased angular momentum, drawn from 
the primary star, which will slow the contraction of the sec- 
ondary Roche lobe. Also, the matter comes from a cool disc, 
so it is likely to come with a reduced specific entropy, en- 
couraging the secondary star to contract. Because of this 
latter effect, and because some mass may be lost to the sys- 
tem rather than re-accreted by the secondary, the period of 
the CV may shorten, following the Roche lobe relation 

P orb oc [± j , (22) 
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where R2 is the radius of the secondary star. However the 
evolution in P orb should be slowed as compared to the stan- 
dard picture. Binaries are therefore expected to spend longer 
at short periods than would otherwise be expected in the 
weak propeller case. 

It would be useful to perform more extensive param- 
eter space explorations and to compare the importance of 
magnetic and accretion torques in different configurations. 
A full population synthesis including binary evolution could 
then be used to asses the precise effect of this weak propeller 
model on the CV period distribution. To achieve the most 
realistic results, an improved model for the magnetic field 
would also be desirable. 
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